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We analyze the onset of classical field configurations after a phase transition. Firstly, we motivate 
the problem by means of a toy model in quantum mechanics. Subsequently, we consider a scalar 
field theory in which the system-field interacts with its environment, represented both by further 
scalar fields and by its own short- wavelength modes. We show that, for very rapid quenches, the 
order parameter can be treated classically by the time that it has achieved its ground state values 
£SJ (spinodal time). 

o 

o : 

(N ! I. INTRODUCTION 

x> : 

■ The emergence of classical behaviour from a quantum system is a problem of interest in many branches of physics 
jjj, from the foundations of quantum mechanics, condensed matter and quantum optics, to quantum computing and 
quantum field theory. In the last few years, experimental evidence supports our theoretical understanding of the 
quantum to classical transition and also opens new avenues in the development of essential tools to understand the 

' , more hidden mysteries of the quantum world Q . 

■ The quantum to classical transition, for a point particle, say, involves two different but very related conditions. The 
first one is that there should be correlations, i.e., the coordinates and momenta of a particle should be correlated in 

f^S • phase space, according to the classical equations of motion. For example, the Wigner function should have a peak 
at the classical trajectories. The second and equally important condition is the elimination of quantum interference 
between these classical trajectories, i.e., decoherence. Once decoherence eliminates interference terms (absent in 
the 'classical world'), the Wigner function becomes a good candidate for a classical probability distribution. The 
generalisation from particles to fields is straightforward, in principle, with the same attributes of correlations and 
' decoherence. 

Oh, The onset of classical behaviour is a natural consequence of a quantum open system, triggered by the interaction 
between the system of interest and its environment. The coupling strength between system and bath sets the timc- 
I scale after which we can consider our system as classical, according to both of the conditions we mentioned before. 
^Slj . This temporal scale is usually called the decoherence time tn- After this decoherence time we do not have macroscopic 
states in coherent superpositions anymore, and a probability distribution can be extracted that evolves by means of 
• 1— 1 . a generalised Fokker- Planck equation. 

Our concern in this paper will be the quantum to classical transition of a single scalar order-parameter field during 
continuous transitions, with the simplest double-well potential. This is an idealisation of the phase transitions that 
- ■ are expected to occur at the GUT and EW scales in the standard big-bang cosmological model || . 

An analysis of phase transitions in quantum field theory that takes the non-equilibrium nature of the dynamics into 
account from first principles has only recently begun to be addressed. In particular, the naive picture of a classical 
order parameter (inflaton or Higgs) field cf> rolling down an adiabatic effective potential, that was once a mainstay 
of cosmological field theory modelling, has been shown to be suspect [Q . Alternatively, the suggestion by Kibble || 
that, while a non-adiabatic approach is essential, causality alone can set saturated bounds on time and distance scales 
during a transition, has been shown to be only partly true. The onset of classical behaviour is absolutely crucial in 
the work of Rcf. || in that the experimental signal for the causal bound is the production of classical defects at a 
prescribed density. Topological defects are inevitable in most transitions, and they may have played a fundamental 
role in the formation of large scale structure (strings) ||. Moreover, superabundance of some topological defects may 
contradict the observational evidence (magnetic monopoles). 



* f.lombardo@ic.ac.uk 

tr.rivers@ic.ac.uk 

tfmazziQdf.uba.ar 



1 



In all the above-mentioned examples there is an order parameter which evolves from the false to the true vacuum 
of the theory: the Higgs fields in GUT and EW phase transitions, the infraton field(s) in inflationary models, etc. 
Although these are quantum scalar fields with vanishing mean value (due to the symmetry of the initial quantum 
state), the order parameter is usually treated as a classical object. Our aim in this paper is to justify this assumption. 

Our approach follows the analysis started by two of us in Ref. §, where we studied the emergence of classical 
inhomogcncities from quantum fluctuations for a self-interacting quantum scalar field. We have investigated there the 
decoherence induced on the long-wavelength field modes by coarse-graining the field modes with wavelength shorter 
than a critical value, in order to show how the system becomes classical due to the interaction with its environment 
(in that case composed of the short- wavelength field modes of the same field). For phase transitions the classicality 
of the order parameter can be analyzed along the same lines by extending the model to accommodate spontaneous 
symmetry breaking. 

This is a difficult problem because, as has been pointed out in the literature, and as we will stress in what follows, 
non-perturbative and non-Gaussian effects are relevant in the analysis of phase transitions. As a trial run for this 
analysis we begin by presenting a toy model in which we will study the spread of a wave packet initially centered 
around the local maximum of a double well potential, paying particular attention to the influence of the environment 
on the Wigner function and on the reduced density matrix. Only later will we will extend our results to quantum 
field theory phase transitions. 

The paper is organized as follows. In the next section we study the evolution of a wave packet initially centered on 
the top of a double well potential. We describe the exact numerical evaluation of the evolution of the wave packet. 
We show that, as the coupling between the system and the environment decreases, the decoherence time increases. 
Due to the nonlinearities of the potential, when the coupling vanishes there is no classical limit, not even classical 
correlations. In section III we analyze the full problem in quantum field theory. Here, the situation is somewhat 
different. Although, again, the decoherence time increases as the coupling to the environment decreases, so does the 
spinodal time, the time for the order parameter to achieve its classical ground-state values. Thus, however weak the 
coupling, decoherence can still occur before the field has fallen to its classical level. 

It is important to say that the issue of how the system evolves into the classical theory has been addressed in Refs. 
|^|,^). For some models, it has been shown that classicality emerges as a consequence of profuse particle creation. A 
non-perturbative large occupation number of long- wavelength modes produces, on average, a diagonal density matrix. 
This dephasing effect occurs at late times. Here we will consider a different model in which classicality is an early 
time event. 

Finally, in section IV, we present our final remarks and conclusions. 



II. TOY MODEL: DOUBLE WELL POTENTIAL AND ENVIRONMENT 



The simplest field theory that permits a phase transition is that of a single real scalar field (/>, with action 



S[<f>] = 



1 



X 

if 



(i) 



with (jj, 2 > 0) Z 2 symmetry breaking. At high temperature the symmetry is restored. As we shall see later, a phase 
transition induced by a sudden temperature quench can be described by an effective field theory in which there is a 
change of sign in the mass term of the scalar field 



(2) 



where m 2 (t) = M 2 > for t < and m 2 (t) = — /i 2 for t sufficiently positive. This change of sign in m 2 (t) breaks the 
global Zi symmetry for positive t. 

Understanding this transition, for even such a simple system, is difficult. As a preliminary exercise, we start by 
considering a toy model composed of a particle, a quantum anharmonic oscillator (the 'system'), linearly coupled to 
an environment composed of an infinite set of harmonic oscillators. The total classical action is given by 

S[x, q n ] = S syst [x] + S cnv [q n ] + Si nt [x, q n ] 

Imh 2 -^ 2 --^) 
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where x and q n are the coordinates of the particle and the oscillators, respectively. The quantum anharmonic oscillator 
is coupled linearly to each oscillator in the bath with strength C n . In analogy with m 2 (f) of ([j]) we consider the simplest 
possible case of an instantaneous quench, in which SIq (t) = f2§ when t < and Oo(f) = — when t > 0. The unstable 
particle, coordinate x, has an initial thermal distribution at temperature T, for which (x) — 0. For t > it finds itself 
on the unstable hump of the potential and falls towards its stable minima at (x 2 ) — 20q/A. 



A. The environment 



In the absence of an environment, and considering A = 0, this situation has been discussed in detail by Guth and 
Pi JTc| , [TT[ ] . In the presence of an environment the relevant objects for analyzing the quantum to classical transition in 
this model are the reduced density matrix, and the associated Wigner function 

p r (x,x',t) = J dq n p(x,q n ,x',q n ,t), (4) 

and 

W x (x,p,t) = ± f dye*" prix + ^x-^t). (5) 



The reduced density matrix satisfies a closed master equation. This has been evaluated by Hu-Paz-Zhang |12| for the 
quantum Brownian motion problem with On positive. Following the same procedure, we can write a corresponding 
master equation for the unstable particle |13| by replacing ilo by ifio in the Hu-Paz-Zhang result 

id t p T (x,x',t) = (x\[H,p T ]\x') 

- ij(t)(x - x')(d x - d x ')p T (x,x',t) 
+ f(t)(x- x')(d x +d x ,)p T (x,x',t) 

— iD(t)(x — x') 2 p T (x, x, t). (6) 

In (||), H = i/syst — \MVL 2 (t) where Cl 2 (t) renormalizes the natural frequency of the particle, ^(t) is the dissipation 
coefficient, and D{t) and f(t) are the diffusion coefficients, which produce the decoherence effects. They depend on 
the properties of the environment as 



n 2 (t) = [ dt' cosh(il t')r){t') 



D(t) = [ dt' cosh(fioi'M*') ( 7 ) 
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where r](t) and v(t) are the dissipation and noise kernels, respectively, 

/■oc 

77(f) = / dwl(u>) sinwf 
Jo 

u(t) = / dujl{uj) coth -— coscjf, 
Jo 2 

and I(u>) — 0(C 2 ) is the spectral density of the environment. 

The first term on the RHS of Eq.(|) gives the usual Liouville-like evolution; the term proportional to 7 produces dis- 
sipation (7 is the relaxation coefficient). The term proportional to the diffusion coefficient D(t), which is proportional 
to [x — x') 2 and positive definite, gives the main contribution to the decoherence since it produces a diagonalization 
of the reduced density matrix. Let us write the reduced density matrix as 

p r \x, x'; f] = G[x, x', f] exp — [x — x 1 ) 2 I D{s) ds . (8) 
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Inserting this expression into the master equation it is easy to to see that the differential equation for G[x,x',t] 
contains the usual Liouville-term plus additional contributions proportional to D, 7, and /. However, none of these 
additional terms is imaginary with the right sign for diffusion. An approximate solution of Eq. (0) is therefore H] 



p T [x,x';t] f» [x,x',t] exp - (x - x') 2 \ D(s) ds , (9) 



where p^ takes into account the unitary evolution. 

Alternatively, one can derive the following evolution equation for the reduced Wigner function of the system |l5f| : 

W T (x,p,t) = {H syst ,W r } PB - \xdl pp W, + 2^{t)d p (pW r ) 

+ D(t)d 2 pp W r - f(t)d 2 px W r . (10) 

Let us concentrate on the evolution equation (|i"o|). The first term on the right-hand side of Eq.([l(]) is the Poisson 
bracket, corresponding to the usual classical evolution. The second term includes the quantum correction (we have 
set h — 1). The last three terms describe dissipation and diffusion effects due to coupling to the environment. In 
order to simplify the problem, we consider a high-temperature ohmic (/(w) ~ lo) environment. In this approximation 
the coefficients in Eq.(|To|) become constants: ^(t) = 70, / ~ and D = 2^^^. The normal diffusion coefficient 
D is the term responsible for decoherence effects and at high temperatures is much larger than 70 and /. Therefore, 
in Eq.(|l(i|), we may neglect the dissipation and the anomalous diffusion terms in comparison to the normal diffusion. 
As we saw in (JoJ) , it is the diffusion exponential in this equation that enforces the approximate diagonalisation of p is 
this coordinate basis. 

It is important to note that this high-temperature approximation is well defined after a time-scale of the order of 
\/{k-QT) ~ 70/-D (with h — 1). The relevant regime of evolution for our systems takes place at times comfortably 
larger than this time-scale. 



B. Numerical Analysis 

We have solved equation ( |io| ) numerically, in the high temperature limit, for different values of the diffusion 
coefficient D, in order to illustrate its relevance in the quantum to classical transition. Details are given elsewhere 
fl6|| . We have chosen as initial condition a Gaussian state centered at xq = Po = with minimal uncertainty (a 2 = 0.5 
and a 2 = 0.5). The Wigner function is initially positive definite, and different from zero only near the top of the 
potential. We have set the coupling constant A = 0.01, the renormalized frequencies Q = u) n = 1 (we are measuring 
time in units of Q) and the bare masses also equal to one. 

It is illustrative to examine first the exact result when the environment is absent (for this case we have solved 
numerically the Schrodinger equation). The initially Gaussian Wigner function begins to squeeze in the x = p 
direction and, before the spinodal time (t sp ~ 2.3) it becomes a non-positive function (Fig. 1). Our definition of 
the spinodal time is that time at which (x 2 )t = 2f2g/A, its symmetry-breaking value. During the evolution, the 
Wigner function covers all the phase space (Fig. 2) and it is clear that it is not possible to consider it as a classical 
probability distribution. Although we started with a special initial state (Gaussian with minimum uncertainty), the 
non-linearities of the potential make the Wigner function a non-positive distribution. 

Let us now consider a coupling with an environment such that the normal diffusion coefficient is D = 0.01. As 
expected, the evolution of the Wigner function is similar to the previous one at early times (Figs. 3 and 4). However, 
as can be seen from Figs. 5 and 6, at long times it becomes positive definite and peaked around the classical phase 
space. 

The effect of the environment is more dramatic for larger values of the diffusion coefficient. In our last example, 
D = 1, the quantum to classical transition takes place almost instantaneously, even before the quantum particle pass 
through the spinodal point (Figs. 7 and 8). 

It is interesting to note that, as the diffusion coefficient grows, the amplitude of the Wigner function decreases. This 
is due to the fact that the decoherence increases with D. The reduced density matrix diagonalizes. As a consequence, 
its 'Fourier transform', the reduced Wigner function, spreads out. 

Our numerical results show explicitly that the existence of the environment is crucial in the quantum to classical 
transition. The decoherence time depends on the temperature and the coupling between system and environment 
through the diffusion coefficient D |M . 
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III. PHASE TRANSITIONS IN FIELD THEORY 



We can now tackle the full quantum 0- field theory of (Q). A simplified analysis of the model has been given elsewhere 
p8| . As before, the onset of classical behaviour is due to the environment. For the infinite degree of freedom quantum 
field, the field ordering after the transition begins is due to the growth in amplitude of its unstable long-wavelength 
modes. For these modes the environment consists of the short-wavelength modes of the field, together with all the 
other fields Xa with which the <j> inescapably interacts (tJ. The inclusion of explicit environment fields Xa is both a 
reflection of the fact that a scalar field in isolation is physically unrealistic, as well as providing us with a systematic 
approximation scheme. 

Although the system field can never avoid the decohering environment of its short-wavelength modes, to demonstrate 
the effect of an environment we first consider the case in which it is taken to be composed only of the fields x a . The 
short-wavelength modes of the 4> field will be considered afterwards. Specifically, we take the simplest scalar classical 
action 

S[<p, X ] = S syat [<t>] + S env [ X ] + Si„tfo x], (11) 
5 syst [0] = j d 4 x |^<W + i/iV - ^0 4 | , 

ScnvfXa] = E / d A X | "SpXa^Xa - T^aXal , 
N „ 

SuA<P,x] = - Ef / d 4 x^(x)xl(x), (12) 

a=l 

where /i 2 , m 2 > 0. 

We see that the quantum mechanics action of (^) is, in large part, a simplified version of© To make any progress 
it is important that the environment be as simple as possible. To that end we have included no explicit self-interaction 
of the Xa fields, which interact through the (f> field as intermediary. However, there is one significant difference in 
that the interactions of the ^-field with the environment in ( [l2| ) are quadratic, and not linear as in ( 3|). The more 
conventional choice of a linear coupling to the environment was made in the previous section (and [16j) to give a 
model directly comparable to similar particle models with no symmetry breaking, for which much work has been 
done. Although it has been adopted for quantum field theories [|ll| a linear-linear coupling is inappropriate and we 
have taken the simplest quadratic-quadratic couplings in S- m t for convenience. The extension to Yukawa couplings is 
straightforward, and will be given elsewhere. We are reminded that, for non-linear couplings like x n qf l in the quantum 
mechanical problem, one expects the master equation to contain terms of the form iD^ n,m \t)(x n — x m ) 2 p r . We shall 
find a similar effect here. The non-linear coupling to the environment is crucial to our conclusions. 

Further, in our present model, the environment has a strong impact upon the system-field, but not vice-versa, 
whenever possible. The simplest way to implement this is to take a large number N 3> 1 of x a fields with comparable 
masses m & ~ fi weakly coupled to the (f>, with A; g a <C 1. Thus, at any step, there are iV weakly coupled environmental 
fields influencing the system field, but only one system field to back-react upon the explicit environment. 



A. Initial conditions 



Before we examine the model in detail, there are some general observations to be made about initial conditions, 
and the way in which the transition is implemented. Like any simple scalar theory the model displays a continuous 
transition at a temperature T c , 

T 2 C = » f 

The environmental fields Xa reduce T c and, in order that T c 2 3> [J. 2 , we must take 

For order of magnitude estimations it is sufficient to take identical g a = g/yN, whereby 
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1 > 1/VN> g ~ A. 



i.e. the Xa 'tadpole' diagrams completely overwhelm the <j) self-interaction tadpole diagram in generating the <f> thermal 
mass. 

With r\ = \J 6/i 2 /A determining the position of the minima of the potential and the final value of the order parameter, 
this choice of coupling and environments gives the hierarchy of scales 

important in establishing a reliable approximation scheme. Further, with this choice the dominant hard loop contri- 
bution of the 0-field to the Xa thermal masses is 

6m 2 T = 0{gT 2 /V~N) = 0(p 2 /N) < p 2 . 

Similarly, the two-loop (setting sun) diagram which is the first to contribute to the discontinuity of the x-field 
propagator is of magnitude 

g 2 T 2 /N = 0(gp?/N 3 / 2 ) « 5m%, 

in turn. That is, the effect of the thermal bath on the propagation of the environmental x-fields is ignorable. 

This was our intention in model-making; to construct an environment that reacted on the system field, but was 
not reacted upon by it to any significant extent. We stress that this is not a Hartree or large-N approximation of the 
type that, to date, has been the standard way to proceed [§|,[l9| for a closed system. 

We shall assume that the initial states of the system and environment are both thermal, at a temperature Tq > T c . 
We then imagine a change in the global environment (e.g. expansion in the early universe) that can be characterised 
by a change in temperature from To to Tf < T c . That is, we do not attribute the transition to the effects of the 
environment-fields. On incorporating the hard thermal loop tadpole diagrams of the \ (and 4>) fields in the <fi mass 
term leads to the effective action for <fr quasiparticles, 

where m^(Xb) = — p 2 (l — Tq/T 2 ) = M 2 > 0. As a result, we can take an initial factorised density matrix at 
temperature To of the form p[To] — p^[To]p x [To], where p^po] is determined by the quadratic part of S°y St [4>] and 
p x [To] by 5env[Xa]- Yet again, the many Xa fields have a large effect on </>, but the 0- field has negligible effect on the 
Xa- Provided the change in temperature is not too slow the exponential instabilities of the 0-field grow so fast that 
the field has populated the degenerate vacua well before the temperature has dropped to zero. Since the temperature 
T c has no particular significance for the environment field, for these early times we can keep the temperature of the 
environment fixed at T x = Tq = 0(T C ) (our calculations are only at the level of orders of magnitude). 

Since it is the system- field <f> field whose behaviour changes dramatically on taking T^ through T c , in this paper we 
adopt an instantaneous quench for T^ from To to Tf — at time t = 0, in which m 2 AT) changes sign and magnitude 
instantly, concluding with the value rn^it) = —p 2 ,t > that we cited in ([l]). An instantaneous quench is sufficient 
to demonstrate the rapidity with which the environment forces the system field to become classical. Meanwhile, for 
simplicity the Xa masses are fixed at the common value m ~ p. 



B. Tracing out the \ fields 

At time t > the reduced density matrix p r [<fi + , <t>~ ■, t] — (<fi + \p r (t)\(f)~) is now 

Pr[</> + >~,£] = / r T>Xa p[(p + ,Xa,<j>~ ,Xa,t], 

where p[4> + , xt i _ iXa : A = ( ( / ,+ Xa |p(£)|0 ~Xa ) i s the full density matrix. Since we would like to be able to distinguish 
between different classical system-field configurations evolving after the transition, we will only interested in the field- 
configuration basis for this reduced density matrix (in analogy with the quantum Brownian motion model of the 
previous section). The environment will have had the effect of making the system essentially classical once p r (t) is, 
effectively, diagonal. [Our earlier comments on dephasing remain valid, and we stress again that our understanding of 



G 



what constitutes classical behaviour is essentially different from the dephasing effects found in Refs. |gU£[]).] Quantum 
interference can then be ignored and the system is said to have decohered. At the same time, we obtain a classical 
probability distribution from the diagonal part of p T (t), or equivalently, by means of the reduced Wigner functional. 
For weak coupling there will be no 'recoherence' at later times in which the sense of classical probability will be lost. 
Its temporal evolution is given by 



Pv[faf , fa ,t] = j d(j)f J dfa J,[fa + ,fa ,t\(j)f,fa ,to] prty+fa ,*o], 

where J r is the reduced evolution operator 

j T [fa+, 4>TMt^<KM = r V4> + I** e^^-^F^^-]. (is) 



The Feynman- Vernon [|TJ influence functional F[cj) + ,<f> ] is defined as 

F[4> + ,<l>~] = J dxti J dx~i PxbctvXliM J rfXaf 



x J™ Vxt J™ V X l exp (i{S cav {xi] + S iat [<t> + , Xt] 

x exp (-j{S" cnv [Xa ] + Sint[<t>~,Xl]})- 



Beginning from this initial distribution, peaked around <p — 0, we follow the evolution of the system under the 
influence of the environment fields, with Hamiltonian determined from (p"l|). From the influence functional we define 
the influence action 5A[(f> + , cf>~} as 

F [(/>+,(/>"] = exp i8A[(f) + , </)-]. (14) 
The total coarse-grained effective action is 

A[cp+, (j)-} = S[<j) + ] - S[<j)-} + 6A[(j) + , 



We will calculate the influence action to one loop (two vertices) for N large using closed-time path correlators. It is 



the imaginary part which contains the information about the onset of classical behaviour. If A = \{<t> +2 — <f> 2 ), we 



find 

ImSA = "fg / d 4 x J d A y A(x)N q (x, y)A(y). (15) 

In ( p^| ) N q (x — y) = KeG 2 ^ + (x,y) is the noise (diffusion) kernel, where G++ is the relevant closed-timepath correlator 
of the X- field at temperature To. Non- leading one-loop terms are smaller by a factor 0(N^ 1 ^ 2 ). 

The first step in the evaluation of the master equation is the calculation of the density matrix propagator J r from 
Eq. (|lg|). In order to estimate the functional integration which defines the reduced propagator, we perform a saddle 
point approximation 



where (jy^ is the solution of the equation of motion S f^+ L \ ^ > + =( j > - = with boundary conditions </>^j(to) = <j>t an< ^ 
4>^(t) — (pf. It is very difficult to solve this equation analytically. For simplicity, we assume that the system-field 
contains only one Fourier mode with k = kg. We are motivated in this by the observation that the long- wavelength 
modes start growing exponentially as soon as the quench is performed and rapidly dominate the fluctuation power 
spectrum Q. Modes with |fco| 2 > p 2 will oscillate. 

Further, we are interested primarily in the order-parameter 

cj)(t) = 0(fc o = 0, t) = lim — / d 3 x(t)(x,t), 

V^oo V J xeV 

and all our subsequent calculations will be for theis zero-frequency mode. We write the spatially-constant classical 
solution as </>(s) = f(s,t) where f(s,t) satisfies the boundary conditions f(0,t) = <j)\ and f(t, t) — fa. Qualitatively, 
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f(s,t) grows exponentially with s for t < t sp , and oscillates for t sp < s < t when t > t sp . We shall therefore 
approximate its time dependence for t < t sp as 

f(s,t) = <piUi(s,t) + (f>{U 2 (s,t), (16) 

where ui(0,i) = 1, ui(t,t) = and «2(0,t) = 0, U2(t, t) = 1, with solution 

s'mh\u(t — s)] . . sinh(us) 

«iM = -A A «a M = ■ r , ■ 

smh(/ii) smh(/ii) 

We shall justify the use of the linear equation later. In fact, we shall see then that it is not an unreasonable 
approximation until almost t sp . 

In order to solve the master equation we must compute the final time derivative of the propagator </ r , and after 
that eliminate the dependence on the initial field configurations <j>^ coming from the classical solutions </rjj. This is 
the usual procedure, see Ref. 0. 

As we are solely interested in the onset of classical behaviour, it is sufficient to calculate the correction to the 
usual unitary evolution coming from the noise kernel. For clarity we drop the suffix f on the final state fields. If 
A = (4> +2 — <fi~ 2 )/2 for the final field configurations, then the master equation for p r {4> + , 4>~ , t) is 

ip r = ^+\{H, Pl }^-)~i^VA 2 D(t)p r + ... (17) 

The presence of A in (|l7j), rather than — <j>_, is a consequence of the quadratic coupling to the environment in S[ n t- 
Since our main concern is with the diagonalisation of p r , it is not necessary to consider the evolution equation for the 
Wigner functional, as has been shown in the literature of quantum Brownian motion [ |H||T6| . The volume factor V 
that appears in the master equation is due to the fact we are considering a density matrix which is a functional of two 
different field configurations, 4>^(x) = ^ , which are spread over all space. The time dependent diffusion coefficient 
D x (t) due to each of the many external environmental \ fields is given by 

D x (t) = 3 [ ds u(s,t) ReG^ + (0;t- s), (18) 
Jo 



where 



u(s, t) 



, ,, u 2 (t,t) 
u 2 (s,t) - — — ui(s,t) 

Ul(t, t) 



(19) 



For the case in hand of an instantaneous quench, u(s, t) = cosh p{t — s) when t < t sp , and is an oscillatory function 
of time when t > t sp . 

Although G++ is oscillatory at all times, the exponential growth of u(t) enforces a similar behaviour on D x (t). In 
the high temperature limit (k^T fi), the explicit expression for that contribution to the diffusion coefficient due to 
the Xa fields alone is 

D x (t) = 5 -^T~ f ds u(s, t) ^ dp 4 cospM, (20) 

where uj 2 = p 2 + m 2 and we have set m 2 = y? . 

For times [it ^> 1, the integration in ( po|) is dominated by the behaviour at s = 0: 

D x (t) ~ (fc B T //i) 2 u(0, t) ~ (fc B T //i) 2 exp[2 M t]. (21) 

The spinodal time t sp is again defined as the time for which (<fi 2 )t ~ if — 6/i 2 /A. For t > t sp the diffusion coefficient 
stops growing, and oscillates around D{t = t sp ). 



C. Short wavelength modes 

In our present model the environment fields Xa are not the only decohering agents. The environment is also 
constituted by the short-wavelength modes of the self-interacting field <fi. Therefore, we split the field as <j> = <^>< + <^>> , 
where the system-field 4> < contains the modes with wavelengths longer than the critical value /U , while the bath or 
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environment-field contains wavelengths shorter than pT l . This gives an additional one-loop contribution D^(t) to the 
diffusion function with the same u(s, t) but a G++ constructed form the short-wavelength modes of the (/>-field as it 
evolves from the top of the potential hill. Without the additional powers of TV -1 to order contributions summation 
of loop diagrams is essential to get a reliable G++. However, it is not necessary to calculate D,p(t) in order to get 
a good estimate of tjj. Since the contribution of D$(t) to the overall diffusion function is positive we can derive an 
upper bound on the decoherence time to from the reliable diffusion functions D x alone. 

However, we would not expect the inclusion of the 0- field to give a qualitative change. Specifically, we note that, if we 
take 0-field propagators dressed by only the simplest tadpole diagrams, and ignore two-loop self-interaction diagrams 
then the diffusion correction due to the </> loop is now similar to that of the x loops. The effect is that the short- 
wavelength modes in the one-loop diagrams from which they are calculated have been kept at the initial temperature 
Tp, on the grounds that passing through the transition quickly has no effect on them. The quench mimics a slower 
evolution of temperature in which only the long-wavelength modes show instabilities that the transition induces. 
That is, with g ~ A and no 1/N factor, the short-wavelength </> modes would have the same effect on the dissipation, 
qualitatively, as all the explicit environmental fields put together. However, at an order of magnitude level there is 
no change, since the effect is to replace g 2 by g 2 + 0(A 2 ) = 0{g 2 ). 



D. The decoherence time 



Using the positivity of we estimate the decoherence time £d for the model by considering the approximate 
solution to the master equation (|l7|), 



Pr [</>< , 0< ; t] < [cj>X ,<£<;*] exp 



VT f ds D x (s) 
Jo 



where is the solution of the unitary part of the master equation (i.e. without environment). It is obvious from this 
(and also from ([l?])), that the diagonal density matrix just evolves like the unitary matrix (the environment has almost 
no effect on the diagonal part of p T ). In terms of the dimensionless fields (f> — (</>< + 0<)/2/i, and 8 — (<f>< ~ (f><)/2p, 
we have T = (1/16)<?V^ 2 <5 2 . 

The system behaves classically when p r is appropriately diagonal. We therefore look at the ratio 



Pr[(f) + S,(f) - 8;t] 



p T [<f>,<f>\t] 



< 



p*[4> + 8,4>-8;t) 



P?[<t>,<P;t] 



exp 



VT I ds D x {s) 

IQ 



(22) 



It is not possible to obtain an analytic expression for the ratio of unitary density matrices that appears in Eq. (|22|) . 
The simplest approximation is to neglect the self-coupling of the system field In this case the unitary density 
matrix remains Gaussian at all times as 



p?$ + 8,$-8;t) 



exph^V 1 ^)], 



(23) 



where p _1 (t), essentially p 2 (4> 2 )^ 1 , decreases exponentially with time to a value 0(A). This approximation can be 
improved by means of a Hartree-like approximation p9[. In this case the ratio is still given by Eq.(23), but now 
p~ 1 (t) decreases more slowly as t approaches t sp . In any case, in the unitary part of the reduced density matrix the 
non-diagonal terms are not suppressed. Therefore, in order to obtain classical behaviour, the relevant part of the 
reduced density matrix is the term proportional to the diffusion coefficient in Eq. ( p2|) , since it is this that enforces its 
diagonalisation. 

The decoherence time tv sets the scale after which we have a classical system-field configuration. According to our 
previous discussion, it can be defined as the solution to 



1 > VT I ds D{s) 
'o 



(24) 



It corresponds to the time after which we are able to distinguish between two different field amplitudes (inside a given 
volume V). 

Suppose we reduce the couplings g, A of the system </>-field to its environment. Since, as a one-loop construct, 
r oc g 2 , X 2 our first guess would be that as g, A decrease, then to increases and the system takes longer to become 
classical. This is not really the case. The reason is twofold. Firstly, there is the effect that T a Tq 2 , and TJ 2 oc A -1 







is non-perturbatively large for a phase transition. Secondly, because of the non-linear coupling to the environment, 
obligatory for quantum field theory, T oc (jr. The completion of the transition finds <f> ~ rf 2 oc A -1 also non- 
perturbatively large. This suggests that T, and hence to, can be independent of A. The situation would be different 
for a linear coupling to the environment, for which (jr would not be present, or a cold initial state in which <f> is peaked 
about <j) = 0. In fact, the situation is a little more complicated, but the corollary that to does not exceed t sp as the 
couplings become weaker remains true. 

In order to quantify the decoherence time we have to fix the values of V, S, and 4>. V is understood as the minimal 
volume inside which we do not accept coherent superpositions of macroscopically distinguishable states for the field. 
Thus, our choice is that this volume factor is 0(fj,~ 3 ) since /z -1 (the Compton wavelength) is the smallest scale at 
which we need to look. Inside this volume, we do not discriminate between field amplitudes which differ by O(fi), 
and therefore we take S ~ O(l). For <j> we set (j? ~ 0(a/X), where A < a < 1 is to be determined self-consistently. 

Note that the diagonalisation of p t occurs in time as an exponential of an exponential. As a result, decoherence 
occurs extremely quickly, but not so quickly that [it <C 1. Consequently, in order to evaluate the decoherence time in 
our model, we have to use Eq. ( pl| ) . On taking the equality in ( pi} ) we find 



g*a a 
whereby 

11 f] 
Ht D — - In TV - - In a ~ ln( ; 



exp[2/^] « * = 0( — ), (25) 



4 2 v T n ^' 



This is an upper bound on tp, but probably still qualitatively correct. The value of a is determined as a ~ y/Ji/Tc 
from the condition that, at time tp, (\<t l \' 2 }t ~ arj 2 . Since a«l, in principle, the field has not diffused far from the 
top of the hill before it is behaving classically. 

For comparison, we find t sp , for which (</> 2 )t ~ ?y 2 , given by 

2 

exp[2/ii sp ] « 0{\). (26) 

The exponential factor, as always, arises from the growth of the unstable long-wavelength modes. The factor T^T 1 
comes from the coth(/3w/2) factor that encodes the initial Boltzmann distribution at temperature T > T c . As a result, 



whereby 1 < fitjj < [xt sp , with 



(27) 



fit sp - [it D ~ lln(-) > 1, (28) 
4 n 



for weak enough coupling, or high enough initial temperatures (we have taken To ~ T c throughout). This is our 
main result, that for the physically relevant modes (with small ko) classical behaviour has been established before the 
spinodal time, when the ground states have became populated. 



E. Back-reaction 



We can now justify our earlier assumption that, for an instantaneous quench, nonlinear behaviour only becomes 
important just before the spinodal time |22]. To see this, we adopt the Hartree approximation, in which the equations 
of motion are linearised so that the field modes ft' now satisfy the equation 



d 2 l2 



M 2 (i) 



fk(t) = o, 



where [n9[ 



d 3 p 



W+4A J -JL. C (p)[f+(t)f-(t) 



(29) 



10 



and C(p) = (coth/3w/2)/2w, uj 2 = p 2 + m 2 . 

A similar result would be obtained by extending our initial 0(2) theory to an 0(M) theory in the large-M limit 
before taking M = 2. 

Our coarse-graining retains only the unstable modes in the integral, which suggests the approximate hybrid 
self-consistent equation for p 2 (t), 

p 2 (t) ~ M 2 - C\ { J^ )m exp ^2 jT df , (30) 

(C = 0(1)), which has the exponential growth of the WKB solution, but non-singular behaviour when /x(t) w 0. The 
exact solution to Eq.(^0|) for i < i sp is /Lt(t) = ^tanh^(i sp — t), irrespective of the values of the temperature T and the 
coupling strength. In fact, we anticipated this in ([27]), when we estimated t sp on the assumption that the backreaction 
would only take effect in the final moments. 

That is, the theory only ceases to behave like a free Gaussian theory with upside-down potential at a time £b, 

t sp -t B = 0{ ( i- 1 ). (31) 

It follows that £b > tjj in our ordering of scales T c 3> fi but, in practice T c needs to be at least an order of magnitude 
larger than p, for this to be the case. 

When ( p8| ) is valid, we see that p r becomes diagonal before non-linear terms could be relevant. In this sense, 
classical behaviour has been achieved before quantum effects could destroy the positivity of the Wigner function W T . 



IV. FINAL REMARKS 



We have shown that, in our model, p T becomes diagonal before the spinodal time at which the order parameter field 
has first populated the ground-state values of the theory. Further, it can also be diagonal before non-linear terms are 
relevant. That is, decoherence can be achieved before quantum effects destroy the positivity of the Wigner function 
W r . Really, our to sets the time after which we have a classical probability distribution (positive definite) even for 
times t > t sp . The existence of the environment is crucial in doing this. Of course, for non-Gaussian or delocalised (in 
the field space) initial states, it is clear that W T will be non-positive definite even in the linear regime, and therefore 
tu should be smaller than the one we evaluated here. In the present work, £d is the classicalisation tame-bound. 

This result goes in the direction of justifying the use of classical numerical simulations for the analysis of the 
dynamics of the long-wavelengths modes after the quench. 
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